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We perform a cosmographic analysis using several cosmological observables such as the luminos- 
ity distance moduli, the volume distance, the angular diameter distance and the Hubble parameter. 
These quantities are determined using different data sets: Supernovae type la and Gamma Ray 
Bursts, the Baryonic Acoustic Oscillations, the cosmic microwave background power spectrum and 
the Hubble parameter as measured from surveys of galaxies. This data set allows to put constraints 
on the cosmographic expansion with unprecedented precision. We also present forecasts for the 
coefficients of the kinematic expansion using future but realistic data sets: constraints on the co- 
efficients of the expansions are likely to improve by a factor ten with the upcoming large scale 
structure probes. Finally, we derive the set of the cosmographic parameters for several cosmological 
models (including ACDM) and compare them with our best fit set. While distance measurements 
are unable to discriminate among these models, we show that the inclusion of the Hubble data set 
leads to strong constraints on the lowest order coefficients and in particular it is incompatible with 
ACDM at 3-<j confidence level. We discuss the reliability of this determination and suggest further 
observations which might be of crucial importance for the viability of cosmographic tests in the next 
future. 

PACS numbers: 98.80.-k, 95.30.Sf, 97.60.Bw, 98.70.Rz 

I. INTRODUCTION 

Cosmology is going through a golden age as we can nowadays start reconstructing the expansion history of the 
universe with unprecedented precision. The huge range of data sets spans a wide realm of observations with het- 
erogeneous nature, providing us with a much more accurate tool for investigating the evolution of the universe 
In particular, most of the observations agree on the evidence that the universe is undergoing an era of positively 
accelerated expansion, requiring the existence of a (more or less conservative) source able to produce it. Since a whole 
set of cosmological models explaining the late time acceleration has been investigated, it goes without saying that 
a sensible test to discriminate among different cosmological evolutions is related to a proper interpretation of high 
redshift data. 

Among these proposed models an important role is played by extended/modified theories of gravity, i.e. by models 
in which the late time acceleration is a by-product of some modified gravitational dynamics (see Ref.[2| for a review). 
In this context, it is clear that the development of a "gravitational dynamics independent" reconstruction of the 
expansion history of our universe does play a crucial role. 

Cosmography provides such unbiased test of the cosmological history by assuming just homogeneity and isotropy 
and then use the so obtained Friedman-Lemaitre-Robertson- Walker (FLRW) metric to express the distances 1 of 
the observed objects as power series in a suitable redshift parameter. The coefficients of such powers, casted into a 
combination of successive weighted derivatives of the scale factor a(t), contain the relevant information for a kinematic 
description of the universe [S-Q ■ 

A comment is necessary here: as already stressed in Rcf.Q the ill-behavior at high z (close and higher than z « 1) 
of the usual redshift expansions strongly affects the results leading in general to an underestimate of the errors. In 
order to avoid these problems, as well as to control properly the approximation associated with the truncation of 
the expansion, it is useful to recast all the involved quantities as functions of an improved parameter y = z/(\ + z) 
[^. HolflH . In such a way, being z £ (0, oo) mapped into y £ (0, 1), it becomes possible to retrieve improved convergence 



Note that depending on which physical quantity one is measuring, it could be more convenient to extract from some data set a particular 
distance indicator than another one. These different quantities have different expressions of the Taylor expansion in redshift, such that 
it could be more natural to estimate cosmographic parameters, whose expression instead does not depend on the analytic expansion, 
in one of these particular frameworks. This ambiguity led also to a misconception about the appropriate definition of distance one 
should investigate. From now on we will refer only to: luminosity distance as the most direct choice in the case of measures of distance 
for Supernovae Type la (SNela) and Gamma Ray Bursts (GRBs); volume distance for Baryon Acoustic Oscillations (BAOs); angular 
diameter distance for the Cosmic Microwave Background (CMB) (see below for their definitions). For a critical discussion about such 
difficulties, see Ref.[|. 
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properties of the Taylor series at high redshift @, UH • 

Some recent papers handle the problem of interpreting the data under a cosmographic perspective using different 
probes [IJ, [Rfl4f5l ] . In this paper we arc going to explore the whole ensemble of data sets and use it to constrain the 
parameters appearing in the expansions of the characteristic scales associated to these indicators: Supernovae and 
Gamma Ray Bursts, Baryon Acoustic Oscillations, Cosmic Microwave Background and Hubble parameter estimates 
(Hub). It is quite obvious that by adding higher order powers to the redshift expansions of such scales it is possible 
to improve the data fitting, since more free parameters are involved. However, for a given data set, there will be 
an upper bound on the order which is statistically significant in the data analysis (Tlj . In this sense it is crucial to 
always determine the order of the expansion which maximizes the statistical significance of the fit for a given data 
set or an ensemble of them. We shall hence determine such order by performing suitable F-tests depending on the 
collection of data sets we shall consider. However, we will show that these statistical arguments can be misleading 
whenever is the case of a model, favored by the -F-test, but with error intervals so large to make the last determined 
parameter completely unconstrained. The parameters obtained within the cosmographic analysis will then finally 
allow to evaluate, in a dynamic independent way, the viability of any theory aiming to explain the current expansion 
of the universe. 

The plan of the paper is as follows: in Section [II] we will describe the cosmographic expansion and the theoretical 
framework (also partly reported in the Appendix), while in Section [IIII wc will present the observational probes used 
for the present analysis. Section [IVI contains the main results using the fourth and fifth order expansions. In Section 
IVl we focus on the forecasts with futuristic data sets while Section I VII is devoted to the study and comparison of the 
cosmographic expansions for different cosmological models. We conclude in Section IVIII 



II. COSMOGRAPHIC EXPANSIONS 



As a pedagogical example, we will discuss first the procedure that has been followed in order to obtain the cosmo- 
graphic expansion for the luminosity distance. As already pointed out, we will start from the only assumption that 
the universe is homogeneous and isotropic, so that the metric describing its properties is the FLRW one 



ds 2 = -c 2 dt 2 + a 2 (t) 



dr 2 



1 — kr 2 



(1) 



using this metric, it is possible to express the luminosity distance g?l as a power expansion in the redshift parameter 
z (or in term of the y-parameter), where the coefficients of the expansion are some functions of the scale factor a(t) 
and its higher order derivatives. 

Following Ref . [lH] . the relation between the apparent luminosity I of an object and its absolute luminosity L defines 
the luminosity distance dh 
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L 
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where r\ is the comoving radius of the light source emitting at time t\ , to is the later time an observer in r = is 
catching the photons, and redshift z is, as usuah defined as I + z = a(to)/a(t\). The radial coordinate r\ in a FLRW 
universe can be written for small distances as 
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with k = —I, 0, + I respectively for hyperspherical, Euclidean or spherical universe. In such a way, it is possible to 
recover the expansion of d^ for small z 



d L (z) = cHo^Z + i(l - q )z 2 - | Tl - qo - 3q 2 + Jo 
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having defined the cosmographic parameters as 
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If instead we use the redshift variable y = z/(l + z) , the definition of the cosmographic parameters will not be affected, 
while now the luminosity distance turns out to be 



d L (y) 




^(<Zo - 3)y 2 + i [11 - 5q + 3q 2 - j + il ko ] y 3 + 0{y 4 ) 



(G) 



where flk = — kc 2 / H^a 2 (to) is the spatial curvature energy density. For a flat universe, flk a = 0. Since we are 
interested in spanning the universe at any redshift, in the following we will use only the formulation of the expansion 
in the variable y. In our analysis we will put constraints up to the cosmographic parameters sq and cq, namely: 
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In the Appendix we list the cosmographic series for the physical quantities involved in our analysis. 



III. OBSERVATIONAL DATA SETS 



In our calculations, we rely on the following current observational data sets: i) SNIa and GRB distance moduli; ii) 
BAO in the galaxy power spectra; iii) CMB; iv) Hubble parameter determinations. 

The SNIa distance moduli provide the luminosity distance as a function of redshift D^(z). In this paper we will use 
the latest SNIa data sets from the Supernova Cosmology Project, "Union2 Compilation" which consists of 557 samples 
and spans the redshift range < z < 1.55 [13] ■ In this data set, they improved the data analysis method by using 
and refining the approach of their previous work [l8j . When comparing with the previous "Union Compilation" , they 
extended the sample with the supernovae from Refs.[l7l Il9j|. The authors also provide the covariance matrix of data 
with and without systematic errors and, in order to be conservative, we include systematic errors in our calculations. 

In addition, we also consider another luminosity distance indicator provided by GRBs, that can potentially be used to 
measure the luminosity distance out to higher redshift than SNIa. GRBs are not standard candles since their isotropic 
equivalent energetics and luminosities span 3 — 4 orders of magnitude. However, similarly to SNIa it has been proposed 
to use correlations between various properties of the prompt emission and also of the afterglow emission to standardize 
GRB energetics (e.g. Ref.[2(|). Recently, several empirical correlations between GRB observables were reported, and 
these findings have triggered intensive studies on the possibility of using GRBs as cosmological "standard" candles. 
However, due to the lack of low-redshift long GRB data to calibrate these relations, in a cosmology-independent way, 
the parameters of the reported correlations are given assuming an input cosmology and obviously depend on the same 
cosmological parameters that we would like to constrain. Thus, applying such relations to constrain cosmological 
parameters leads to biased results. In Ref.[2l[ this "circular problem" is naturally eliminated by marginalizing over 
the free parameters involved in the correlations; in addition, some results show that these correlations do not change 
significantly for a wide range of cosmological parameters [22|, [23| . Therefore, in this paper we use the 69 GRBs over 
a redshift range z <E [0.17,6.60] presented in Ref.(23|, but we keep into account in our statistical analysis the issues 
related to the circular problem that are more extensively discussed in Ref . [2l[ and also the fact that all the correlations 
used to standardize GRBs have scatter and a poorly understood physics. For a more extensive discussion and for a 
full presentation of a GRB Hubble Diagram with the same sample that we used we refer the reader to section 4 of 
Ref. 

In the calculation of the likelihood from SNIa and GRBs, we have marginalized over the absolute magnitude M 
which is a nuisance parameter, as done in Refs.[24l. [25j 
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BAOs have been detected in the current galaxy redshift survey data from the SDSS and the Two-degree Field Galaxy 
Redshift Survey (2dFGRS) f3lT[33l |. They can directly measure not only the angular diameter distance, Da{z), but 
also the expansion rate of the universe, H(z), a powerful tool for studying dark energ y |3 4ll. Since current BAO data 
are not accurate enough for extracting the information of Da{z) and H{z) separately [35j . one can only determine an 
effective "volume" distance [3l| 



D V (Z) 



1/3 



(10) 



In this paper we use the Gaussian priors on the distance ratio of the volume distances as recently extracted from 
the SDSS and 2dFGRS surveys [H[ at z = 0.35 and at z = 0.2 (the two mean rcdshifts of the surveys) 

^p^f = 1.736 ±0.065 (1.). (11) 
Dy(z = 0.2) 

The x 2 °f BAO data used in the Monte Carlo Markov Chain analysis will thus be 

a _ (Dy(z = 0.35)/£>y(z = 0.2) - 1.736) 2 

XBAO " ^2 ■ W 

It is worth stressing here that actually both the physics and the data of BAOs depend on the content in matter of 
the universe f2 m . Hence, they arc a priori dependent on a chosen dynamical framework (see also Ref.[36| for a review). 
This issue is usually ignored in the data analyses performed in the literature. However, such an approximation turns 
out to be valid if one does not range far away from the typical fiducial model firstly used in the determination 
of the physical data points. Indeed, the deviation of different models from the fiducial one can be parametrized 
and estimated by the ratio Z?y(new model) / Dy (fiducial model). The impact of the spacetime priors on the power 
spectrum measurement was analyzed in ref. [37| and led to the conclusion that the ratio Eq. ([lip is only weakly 
dependent on dynamical features. Hence, we can safely use BAOs as a further tool to constrain the cosmographic 
parameters. 

Next step in our analysis is the inclusion of the CMB measurement which is sensitive to the distance from the last 
scattering surface via the locations of peaks and troughs of the acoustic oscillations. This data constrains the curve 
of the cosmological history at very high redshift, z ~ 1100, and hence could be very helpful for discriminating among 
competing theoretical models for dark energy, as they necessarily have to coincide at z < 1 - see for example ref. 
[45j . The sound horizon at decoupling 2 , r s (z#), sets a physical scale for the baryon-photon oscillations depending on 
the baryon density, the photon energy density, and the cold dark matter density. Now, it is known that the angular 
diameter distance Da(z) describes the ratio between the proper size of an object at a certain redshift z and the 
correlated observed angular size. The angle 6a under which the sound horizon is observed today is given by 

6 A = irl^ 1 = r s (z*)/D A (z*) = 0.593° ± 0.001° (1 a) , (13) 

where I a denotes the location of the first peak in the multipole space. As for BAOs, the dependence on the cosmological 
density parameters would not allow the use of CMB observables in a fully cosmographic approach. However, following 
Ref. [38| is possible to give model- independent cosmological constraints if one specifies some extra physical assumptions 
to be fullfillcd by cosmological models. The CMB power spectrum today (apart from the low multipolcs shape) is 
shared by models having the same primordial perturbation spectra and the same value of ^cdm and fibaryon- For 
this reason a basically mo del- independent approach can be pursued by restricting our analysis to models having a 
standard physics up to the decoupling era; asking that new physics after decoupling only modifies the small angle 
spectrum by changing the overall amplitude and Da (2*), and requiring that any multipole-dependent effect at late 
time remain small. Such assumptions, while is cutting away some models like f(R) models with no Dark Matter (39j 
or models with new radiation degrees of freedom, are still general enough to cover most of the cosmological models 
on the market. 



2 In our calculation, we choose z, = 1091.3, the best fit value obtained by the WMAP group |46|| . 
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Finally, we add the direct determinations of the Hubble parameter H(z) to constrain the cosmographic expansion. 
Since the Hubble parameter depends on the differential age of the Universe as a function of redshift, 



measuring the dz/dt could directly estimate H(z). Ref.|2a| used the Sloan Digital Sky Survey data and obtained a 
measurement of H(z) at the redshift z ~ 0. In Ref.[13|, the public data of Gemini Deep Survey (GDDS) survey [28| 
and archival data [29| were used in order to get the differential ages of galaxies. In practice, they selected samples 
of passively evolving galaxies with high-quality spectroscopy, and then used stellar population models to constrain 
the age of the oldest stars in these galaxies (we refer to their paper for a more exhaustive explanation of the method 
used). After that, they computed differential ages at different redshift bins and obtained eight determinations of the 
Hubble parameter H(z) in the redshift range z G [0.1, 1.8]. We calculate the x 2 value of this Hubble parameter data 
by using 



where H th (z) and H ohs (z) are the theoretical and observed values of Hubble parameter, and cfh denotes the error 
bar of observed data. We also make use of the newly released prior on the Hubble parameter Hq, which consists 
of a measurement of the Hubble parameter obtained by the Near Infrared Camera and Multi-Object Spectrometer 
(NICMOS) Camera 2 of the Hubble Space Telescope (HST). 

These observations fix the parameter Hq = lOOft-o (km/s)/Mpc by a Gaussian likelihood function centered around 
i?o = 74.2 (km/s)/Mpc and with a standard deviation a = 3.8 (km/s)/Mpc [3(|. We stress that all the mentioned 
methods for determining H(z) are basically "gravitation theory independent". 

An important point must be underlined: the Taylor series of the Hubble parameter already includes into the 
coefficient of the n-th y-power the same number of cosmographic parameters of the other series expanded up to the 
(n + l)-th y-power (see Appendix). This is due, in comparison with the other distance definitions above, to an extra 
derivative with respect to time included in the definition of the Hubble parameter (see also Ref.[l5j). 

For this reason, and for the different nature of the Hubble data, we will initially consider constraints (based on 
standard candles and rulers) of the expansion coefficients associated to different notions of distances; at the end, we 
will add the Hubble data using one order less in the y-power expansion with respect to the distance data in order to 
constrain the same set of parameters. 

In order to compute the likelihood, we use a Monte Carlo Markov Chain technique as it is usually done in order 
to explore efficiently a multi-dimensional parameter space in a Bayesian framework. For each Monte Carlo Markov 
Chain calculation, we run four independent chains that consist of about 300,000 — 500,000 chain elements each. We 
test the convergence of the chains by using the Gelman and Rubin criterion [Z(| with R — 1 of order 0.01, which is 
more conservative than the often used and recommended value R — 1 < 0.1 for standard cosmological calculations. 



In this section we present our main results on the constraints for the cosmographic expansion from the current 
observational data sets. 

With the accumulations of new data and the improvements of their quality, it is of great interest to estimate the free 
parameters in the polynomial terms of highest order. In our previous paper [TTj we already showed the inconsistency 
of the results in the analysis of the cosmographic expansion caused by early truncations of the power series. For these 
reasons we will now present the results obtained for the most meaningful term of the expansion. In order to find out 
which is the most viable truncation of the scries for a given data set, one can use a test comparing two nested models 
(in this case, two different truncations of the Taylor series). The F-test provides exactly this criterion of comparison, 
identifying which of two alternatives fits better, and in the more statistically significant way, the data. 

In such test, one assumes the correctness of one of the models (the one with less parameters), and then assess the 
probability for the alternative model to fit the data as well. If this probability is high, then no statistical benefit 
comes from the extra degrees of freedom associated to the new model. Hence, the smaller the probability, the more 
significant the data fitting of the second model against the first one will be. Quantitatively, the F-ratio among the 
two polynomials is defined as 




(15) 



IV. DATA ANALYSIS 



F EE 
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(16) 
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TABLE I: Constraints on the cosmography parameters up to fifth order term from different data combinations. 



Data 


SNIa 


Parameter 


<7o 


jo 


so 


Co 




Best Fit 


-0.41 


-1.99 








Mean 


-0.41 ± 0.16 


-1.99 + 1.36 








Xmin/d.O.f. 






549.69/555 






Data 


SNIa+GRB 


Parameter 


qo 


jo 


so 


Co 


H 


Best Fit 


-0.78 


5.03 


50.18 






Mean 


-0.76 ± 0.26 


4.82 + 4.07 


53.57 ± 46.38 






Xmin/d.O.f. 






628.70/623 






Data 


SNIa+GRB+BAO+CMB 4 t/j order 


Parameter 


qo 


jo 


so 


CO 


Ho 


Best Fit 


-0.32 


-2.57 


-18.40 






Mean 


-0.28 ±0.17 


-2.88 ± 1.64 


-17.61 + 2.56 






Xmin/d.O.f. 






633.33 / 625 






Data 


SNIa+GRB+BAO+CMB 5" 1 order 


Parameter 




jo 


so 


CO 


Ho 


Best Fit 


-0.17 


-6.92 


-74.18 


-10.58 




Mean 


-0.49 ± 0.29 


-0.50 ± 4.74 


-9.31 + 42.96 


126.67 + 190.15 




Xmin/d-O.f. 






627.61/624 






Data 


SNIa+GRB+BAO+Hub+CMB (5 th 


order) +Hub (4 th 


order) 


Parameter 


qo 


jo 


so 


CO 


Ho 


Best Fit 


-0.24 


-4.82 


-47.87 


-49.08 


71.65 


Mean 


-0.30 + 0.16 


-4.62 ± 1.74 


-41.05 + 20.90 


-3.50 + 105.37 


71.16 + 3.08 


Xmin/d.O.f. 






639.81/633 







where N is the number of data points, and represent the number of parameters of the i-modcl. The P- value, i.e. 
the area subtended by the ^-distribution curve delimited from the F-ratio point, quantifies the viability of matching 
models as already mentioned. We use the threshold of 5% as the significance level on the P-value under which the 
model with one more parameter fits the data better than the other one. 

We already found in Ref.[ll[ that variations of the total energy density of the universe Qo = 1 — Oft , with the 
spatial curvature parameter ranging between -1 and 1, have a negligible effect on the cosmographic constraints. This 
is basically due to the fact that the error bars for the cosmographic parameters are still quite large in comparison 
with the best fit values. Nonetheless, it is worth noting that this will not be necessarily the case when future data, 
especially at moderate or high rcdshift, will substantially improve the constraints. It is then possible that future 
cosmographic analysis might have to include the spatial curvature effects in the reconstruction of the overall history 
of the universe. This would be the cosmographic expansion counterpart of the strong sensitivity on f2fc showed by 
the reconstruction of w(z) [4lj . 

We then assume f2fc = in our analysis and only present the results for the cosmographic parameters, instead of 
their combinations with r2^ , since the effect of curvature can be for the moment safely neglected. Table U shows the 
constraints on the cosmography parameters as obtained from different data combinations. 

We start performing the data analysis with the SNIa data only. We find that already at the fourth order term 
in the expansion, the minimal x 2 is 549.59. This is not reduced significantly when compared with the constraint of 
the third order case, which has Xmin = 549.69. Hence, introducing the snap free parameter sq does not improve the 
constraints. Indeed, using the F-test, we find a F-ratio of 0.11 and a P-value of 73.93%. Therefore, cosmography up 
to the fourth order term does not fit the SNIa data significantly better: the cosmographic expansion up to the jerk 
term jo (third order) is enough. 

After adding the GRB data, the fourth order case could give a better constraint than third order only. When 
comparing the SNIa+GRB results, the minimal x 2 nas been reduced by about five (x 2 = 628.70 instead of \ 2 = 



FIG. 1: One-dimensional likelihood distributions for the parameters go, jo, so and co for the data combinations 
SNIa+GRB+BAO+CMB+Hub. 



633.32). Using one more time the F-test to contrast third and fourth order expansions, we find F-ratio = 4.59, 
P-value = 3.26%. Thus in this case the fourth order term indeed helps to fit the observed data significantly better. 
The inclusion of GRBs was found to constrain the deceleration parameter q as q = —0.76 ±0.26, so confirming that 



our universe is undergoing an accelerated expansion with a confidence level which is marginally at 3 a [111 ] . The 
-F-test does not suggest to further improve the expansion up to fifth order. 

Including the data point related to BAO does not improve significantly the constraints. The constraining power of 
BAO is rather weak since there is only one BAO data point and its redshift is much lower than those of SN and GRB 
data. For this reason we will consider directly the data set that includes both BAO and CMB. 

When the CMB angle 8a defined in Eq. (fl"3f is added into our analysis, the cosmographic curve is constrained at very 
high redshift, ~ 1100 (namely y ~ 1). Even though the CMB observable is providing just one data point, due to its 
high redshift it is in principle very helpful for discriminating among competing theoretical models producing late time 
accelerated expansion, since these necessarily converge to the same cosmological history at small z (see for example 
Ref.[45[). The large difference between the two Xmin °f the f° ur th and fifth order expansion in powers of y of the 
distances, implies that the latter is the (statistically) more reliable parametrization (Fratio = 5.69 and P V aiuc = 0.02%), 
giving a result very close to the ACDM prediction. Sixth order expansions does not give any substantial statistical 
improvement. 

As already stated at the end of the previous section, Hubble data must be added and analyzed cautiously, since 
they are inhomogencous with respect to the previous data sets both in nature and mathematical handling. In Table 
I we present directly the results for the constraints up to the cq cosmographic parameter, since this truncation turns 
out to be strongly favored with respect to the previous nested model (F ra tio = 19.77 and P V aiue < 0.01%). 



3 Note that our best fit here is different from the one reported in Ref. due to our use of the improved SN catalogue "Union2 
Compilation" . 
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FIG. 2: Theoretical predictions of distance moduli (left panel) and Hubble parameter (right panel) from the best fit model with 
the full data combination, together with the observed data sets. We also show the curves obtained in the ACDM framework 
for comparison (thin black solid lines). 



The theoretical curves of /j,(z) and H(z) are in good agreement with the observed cosmological data, as shown in 
figured The constraint on H is close to the usually quoted value, namely at 68% confidence level is H = 71.16±3.08 
(km/s)/Mpc. One can sec that the addition of the Hubble data leads to relevant improvements in the determination 
of the other cosmographic parameters with the exception of cq, which is still basically unconstrained. We also 
checked whether the next cosmographic parameter had to be included. We find that, for the richest combination 
SN+GRB+BAO+CMB+Hub, the new Xmim i s extremely close to the value in Table HI Therefore, we stop our 
analysis here. 

It is here interesting to underline the power of the y-expanded series (convergent as long as y < 1) allowing us to 
describe the whole cosmological history with the use of relatively few parameters. This circumstance becomes for 
example evident in the left panel of figure O where the furthest GRB data point reaches the distance of y ~ 0.87 
(corresponding to z ~ 6.6). 



V. FORECASTING 



Since the present data do not give yet very stringent constraints on the parameters of cosmography, especially for 
the parameter of fifth order term, it is worthwhile discussing whether future data could determine these parameters 
more effectively For this purpose in what follows we shall perform new analysis of possible future constraints by 
choosing, as a fiducial model, the best fit parameter set for the cosmographic expansion up to the fifth order term as 
fixed by the combination of all the previously considered data sets. 

The projected satellite SNAP (Supernova / Acceleration Probe) would be a space based telescope with a one square 
degree field of view with 10 9 pixels. It aims at increasing the discovery rate for SNIa. to about 2000 per year in the 
redshift range 0.2 < z < 1.7. In this paper we simulate about 2000 SNIa according to the forecast distribution of 
the SNAP. For the error, we follow the Ref.[42| which takes the magnitude dispersion 0.15 and the systematic error 

cr sys = 0.02 x z/1.7. The whole error for each data is given by cr ma ,g(zi) — \J a sy S { z i) + 0.15 2 /rii, where rij is the 

number of Supernovae of the i-th redshift bin. Furthermore, we add as an external data set a mock data set of 400 
GRBs, in the redshift range 0.4 < z < 6.4 with an intrinsic dispersion in the distance modulus of = 0.16 and with 
a redshift distribution very similar to that of figure 1 of Ref . |43| . 

Regarding a future BAO data set, we adopt the predicted performance of the BOSS survey in SDSS III, which will 
measure the angular diameter distance cIa(z) and the Hubble expansion rate H(z) of the Universe over a broad range 
of redshifts. The measurement precision for cLa{z) is 1.0%, 1.1%, and 1.5% at z = 0.35, 0.6, and 2.5, respectively, and 
the forecast precision for the H(z) is 1.8%, 1.7%, and 1.2% at the same redshifts j44[. We also impose a Gaussian 
prior for the current Hubble parameter Hq with the error of 1% provided by a future direct measurement. 

Next coming CMB measurement, mainly via the Planck satellite, could give quite accurate constraints on the 
cosmological parameters. The error bar of 9a could be shrunk by a factor of 3, namely, the standard derivation 
(i0 A = 0.0003°. 

Using all these future mock data, we get the standard derivations of cosmographic parameters: cr qo — 0.02, <jj ~ 
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0.08, fj Sa — 1.95, a Ca = 14.20 and <jhq = 0.48, respectively. We hence see that the constraints on the parameters 
provided by the future mock data can be strongly improved in comparison with the current constraints in Table HI 

VI. COSMOGRAPHIC SELECTION OF VIABLE COSMOLOGICAL MODELS 

In the case of the standard flat ACDM model (namely a model described by Cold Dark Matter with the adding of 
a cosmological parameter) the set of cosmographic parameters results to be (up to fifth order) 



3 
2 

J'o = 1 



10 — n"m — 1 



s — 1 — ^^m 

c = l + 3fi roo + yfi^ o . (17) 

We can use independent probes to constrain the free parameters of the cosmological model, in this case, for example, 
the WMAP estimates of fl mo for the ACDM model. 

The theoretical predictions of the cosmographic parameters in the standard ACDM model are: qo = —0.588, jo = 1, 
so = —0.238 and Co = 2.846, where we set the current matter density to be the best fit value fi mo = 0.275 4 obtained 
by the WMAP group j4f|. Future experiments, in this perspective, will give stricter constraints on the validity of 
such hypothesis. 

Another example is provided by the Dvali-Gabadadze-Porrati (DGP) self-accelerating braneworld model [47[. The 
presence of the infinite- volume extra dimension modifies the Fricdmann equation as: 



^ = n fe (i + zf + (v^T+ y/n rc + n mo (i + z) 3 ) 2 , (is) 



iF 
Ml 

with £7 ?v = l/4r^i?Q accounting for the fractional contribution of the bulk-induced term with respect to the crossover 
radius r c . In a spatially flat universe, = and tt rc = (1 — f2 mo ) /4, the previous equation reads 

(19) 



H 2 



Hi 




+ fi mo (l + z) 3 



so, expanding both the sides of Eq. (|19[) and equating terms of the same power, we obtains the following expressions 
for the cosmographic coefficients as functions of the free parameter Q mg (see also Ref.[48j]) 

-1 + 2£l mo 

qo = 



1 ~h ^mo 



Jo 



i + 3tt„ l0 - m 2 mo + ion 3 mo 



(l + ^mo) 3 

i - m mn + una - i34fii„ + 86r>i„ - sonl 



so 



co = 



(l + ^no) 5 

1 + Un mo - Um 2 no + 1259^- 1996^+ 3828^- 16040^+ 880^, 

(l + O mo ) 7 



(20) 



4 The estimate of O mo is of course known within a certain error. Prom now on, for illustrative purposes, we will retain the best fit values 
of the free parameters, independently estimated in every single cosmological model, as the fiducial ones, without taking into account 
the associated errors. 
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In Ref.[49(, the DGP model has been constrained starting from gravitational lensing statistics; considering the 
fractional amount of matter obtained therein, f2 mo = 0.30, we obtain the following set of values for the previous 
parameters: q Q = -0.308, j = 0.742, s„ = -0.432, c = 2.926. 

We will now take into account the so-called Cardassian cosmology [Hoj], a model whose modification with respect 
to standard ACDM cosmology consists in the introduction of an additional term p n in the matter source of the 
Fricdmann equation, so that now it can be written in term of the fractional matter density as: 

£ = O mo (l + zf + (1 - Q mo ) (1 + zf n . (21) 

Performing one more time the expansion of both sides of the equation, the first four cosmographic parameters can 
now be expressed as functions of the two parameters f2 mo and n 

qo = \ + \{i-n){n mo -i) 

Jo = \ [2 + 9n(fi mo -l) + 9n 2 (l-a™ )] 

so = \ [4 - 180 mo - 9n (4 - 7Q OTo + 3Q roo 2 ) + 9n 2 (ll - 17fi mo + 6SI. 
-27n 3 (3 - 4f2 mo + fi„ i0 2 )] 

- ^l + 30 mo + — — + lg 

-- A (32 - mn mo + 29m 2 mo - 6480^ + 405O£J n + 

+ ^ (136 - 242Q mo + 349ft 2 io - 648^ + 405fi*J n 2 
27 



mo 



._ - 73Q mo + 54^ - 720^ + 45n^J n 3 + 
81 

+ (39 - 560 mo + 26tt 2 iQ - 240^ o + 15fi£J n 4 , (22) 



and for the referring values n mo = 0.271 and n = 0.035 |51|, Eq.((22J) reads g = -0.555, j = 0.890, s = -0.384, 
c = 3.660. 

Finally, we want to show the coefficients in the cosmographic approach of the CPL parametrization [52| for the 
equation of state of Dark Energy. If we suppose to be in a flat universe, then the Friedmann equation is: 

% = n mo (i + zf + (i - n ma ) (i + 2 )3(i+-o+^) e -^ t (23 ) 

and the related cosmographic terms result to be (confront also with ref. [53| ) 

3 

Qo = 1 + ^^o (1 - S^mo) 
3 

Jo = ( 3u, o + 3wq + w a ) (1 - Q mo ) 

7 33 9 

SO = _ "J" (1 ~ fi m ) Wa ~ | (1 - fimo) [9 + (7 - O mo ) lU a ] H>0 - 

9 27 

-- (i - ^ mo ) (i6 - 3f> mo ) w 2 - — (i - n mo ) (3 - n mo ) «,g 

c = -(70 + 3w (-71 + 3u; (-7 + Q mo ))(-l + O roo )) + 
+ ^ (-1 + Q mo ) (-163 + 3w a (-82 + 2m mo )) w + 
+ 9 - (-1 + Q mo ) (-134 - 69w a + 3 (14 + llt« ) TOo ) w 2 + 
+i (1269 - 1917n roo + 648^) ^ + I ( 486 _ 8 i Q mo + 32451^) w 4 ; 

(24) 
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TABLE II: Comparison among cosmographic parameters of different cosmological models. For every model, the evaluation of 
the cosmographic parameters, for a pedagogical issue, is based on the best fit values of the free parameters introduced in the 
dynamics and measured with independent probes. However, the value of the jerk parameter for ACDM model is an exact value, 
see also (17). The values of the cosmographic parameters are compared with our best fits for the series truncations studied 
in the last two lines of Table I. The distances indicators data set includes SNIa, GRB, BAO, CMB. The complete data set is 
obtained adding Hubble data. 



Parameter 


qo 


jo 


so 


CO 


ACDM 


-0.588 


1 


-0.238 


2.846 


DGP 


-0.308 


0.742 


-0.432 


2.926 


Cardassian 


-0.555 


0.890 


-0.384 


3.660 


CPL Parametrization 


-0.511 


0.342 


-2.260 


1.383 


Best fit 

SN+GRB+BAO+CMB (5 th order) 
SN+GRB+BAO+CMB (5 th order) + Hub (4 th order) 


-0.49 ±0.29 
-0.30 ±0.16 


-0.50 ±4.74 
-4.62 ± 1.74 


-9.31 ±42.96 
-41.05 ± 20.90 


126.67 ± 190.15 
-3.50 ± 105.37 



assuming the values suggested by the seventh- year-release of WMAP [46[ for the three free parameters, f2 TOQ = 0.275, 
w = -0.93 and w a = -0.41, wc get q = -0.511, j = 0.342, s = -2.260, c = 1.383. Table HI shows the values of 
the cosmographic parameters in the different models taken into account. 

It is interesting to note a couple of issues in the comparison of cosmological models with our best fits. Firstly, 
the best fit for the cosmographic parameters of the SN±GRB±BAO±CMB data set is perfectly compatible with 
the estimates of the cosmological parameters for a broad variety of models. However, it is easy to realize that the 
currently available data sets would not allow yet to distinguish among the different cosmological models. In fact, Table 
HTl shows that the error bars are still too large with respect to the differences among the cosmographic parameters 
of the cosmological models. Nonetheless, the previously discussed forecasted improvement in the quality and the 
quantity of such data (ameliorating by at least a factor ten the error bars on the cosmographic parameter) should be 
able to discriminate among competing models. 

On the contrary, the best fit of the widest combination of data (that is with the inclusion of the Hubble parameter 
determinations via the differential age technique), seems to exclude, at a 3-cr around the jerk mean value jo, almost 
all cosmological models, including ACDM (with the only exception of the CPL modelization, that is still marginally 
compatible) . We have already discussed the intrinsic difference of the Hubble data and why their use should be taken 
cautiously. It seems clear that this data set while being very powerful in reducing the error bars, is simultaneously 
introducing strong deviations from the mean values determined via standard candles and rulers. This puzzling 
discrepancy in the results does not seem related to the order of the truncation: wc observed a similar behavior even 
for (statistically not favored) early or late truncations of the series. 

However, it is also true that the high rcdshift measurements of the Hubble parameter are based on fitting galaxy 
spectra. As such, this data set is strongly dependent on this fitting procedure which may introduce systematic effects. 
For this reason, we deem estimates based on the Hubble data currently less robust than those based on standard 
candles and rulers. Nonetheless, their use here serves to show the possible key role these data could play in the 
future of Cosmography as they appear to be very effective in reducing the error bars and very sensitive tracers of 
the cosmological history. We hence conclude that our analysis strongly suggests further investigation of this apparent 
tension between the Hubble data and ACDM (and many competing models) via a refinement of the determination 
methods developed in Ref s . [26l [27| . 

VII. DISCUSSION AND CONCLUSION 

Reaching the highest possible redshift allowed by data is a fundamental tool to discriminate among competing 
cosmological models. Given that most of the models are built in order to recover Dark Energy at low redshift, their 
expansion histories are degenerate at late times. To break such a degeneracy, it is required an improvement on the 
knowledge of the early universe expansion curve: this aim can be achieved only by an accurate determination of the 
higher order parameters, and higher terms in the cosmographic expansion can be consistently reached only using 
(very) high redshift data. 

Compared to our previous work we have added Baryonic Acoustic Oscillations (that are distance indicators at 
z ~ 0.3), a compilation of high redshift Hubble parameter measurements and, at least for a wide gamut of cosmological 
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models, CMB data about the angular size of the sound horizon. This improved data set is helpful in that, apart from 
better constraining the previously studied cosmographic parameters, it also allows to cast constraints on the next 
order, so far unbound, expansion coefficient. 

The analysis is performed by using Monte Carlo Markov Chains in the multidimensional parameter space to derive 
the likelihood. As a first step, we consider the most recent catalogs of standard candles, namely Supernovae Type la 
and (properly standardized, see discussion in Section [illj) GRBs. A combination of such data gives constraints up to 
the so parameter in the cosmographic series. We have also used the BAO (albeit the mildly improve the cosmographic 
series fitting) discussing the reliability of such tools in this context. 

Secondly, we add data at higher redshift from different probes to further improve the constraints. The CMB data 
account for a very stable and well determined scale. It is worth noting here, anyway, that on the contrary of the 
other probes, CMB data provide the problem of a lack of universality in the cosmographic approach. Unfortunately, 
the set of parameters extracts from CMB observations is not truly independent from the dynamics of the underlying 
gravitational theory. Its definition, in fact, strictly depends on the assumption of a cosmological model that behaves 
as General Relativity plus a content of matter of arbitrary nature. It is hence impossible to use it straightforwardly 
within a purely cosmographic analysis which wants to apply also to non-standard cosmologies (based on exotic modified 
gravity theories) 5 . In this paper we proposed CMB data constraints on the cosmographic series by restricting the 
results to a slightly smaller variety of models. A desirable full solution to this problem can be achieved "standardizing" 
somehow the CMB parameters or alternatively identifying other CMB observables which could be used as standard 
rulers (at least approximately, as for BAOs). We leave this to future investigations. 

We then added the high redshift measurements of the Hubble parameter. We found that thanks to these data and 
the CMB one, it is possible to ameliorate the knowledge of the cosmographic expansion up to the Co parameter. 

As a completion of our analysis, we have also discussed foreseeable constraints from futuristic data sets provided 
by projected experiments. We showed that a strong reduction of the typical errors on the parameters estimates is a 
realistic goal: future surveys, indeed, do have a solid possibility to sufficiently reduce the uncertainties on the lowest 
order parameters by a factor ten at least, gaining a concrete chance to assess the viability of alternative cosmological 
models (possibly based on different dynamics). 

Finally, we calculated the cosmographic parameter sets for a sample of cosmological models with alternative dy- 
namics (using the so far available best fits for their free parameters). We showed that, while the data set including 
"standard" distance indicators gives a best fit with which all the cosmological models are still compatible, the in- 
clusion of the Hubble data introduces a tension between the observed cosmographic parameters and the parameters 
calculated for different models and in particular with ACDM which appears to be ruled out at 3-er due to the jerk best 
fit value. We have discussed the reliability of such observation taking into account the inhomogeneity of the Hubble 
data set with respect to the distance indicators ones. While there might be systematic uncertainties in this data due 
to their complex determination method, we stressed that our analysis strongly suggests they might play a key role in 
reducing the errors on the estimates of the cosmographic parameters and hence in making Cosmography effective in 
discriminating among competing cosmological models and gravitational theories. 

In conclusion, the search for high redshift standard rulers and most of all the improvement of the data coming 
from galaxy surveys seem to be what could possibly bring cosmographic studies into a mature stage and make them 
powerful, gravity theory independent, tools for selecting among theoretical scenarios. We hence hope that this work 
will further strengthen the case for proposed experiments aimed at improving our knowledge of the cosmic evolution 
of the high redshift universe (e.g. Ref.j55|). 



Appendix: Series Expansions 

We present here more extensively the expansions used in this work. A flat universe, k — 0, is assumed in all the 
expressions below. 
Hubble parameter: 



H[z(y)} = H 



l + (q + l)y + y 2 



2 30 % 



+ Qo + 1 + 



5 Of course, CMB observables can be used within a given gravitational dynamics to fix the free variables of a cosmological model [54| 
and hence calculate the corresponding cosmographic parameters to be confronted with those determined purely on the base of standard 
candles and rulers, as wc also showed as application to the evaluation of cosmographic parameters in several cosmological models. 
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+ gy 3 (-4j g + 3j + 3gg - 3gg + 6g - s + 6) + 
+ 24 y4 ^ _4j ° + 25j09 " ~ 16joq ° + 12jQ + C ° ~ 15<7 ° + 



Luminosity distance: 
dL[z(y)} 



1 



3 
2 

5 jo go 
12 

ii_ 

12 



go 
2 



24 

7j g ; 



5?i 



3 jo 9o 



1 On 2 


f 7goSo 


+ 24g - 4s 4 


6~ 


JO 

~ 6~ 


11\ 


7g 2 


13go 


S+§) + 


8 


12 


47j 




7g 1 9g 3 


120 


— — H 
120 


8 8 




goSo 


77g 3s 




8 


60 40 



0[y 5 }; 



(25) 
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Angular distance: 

d A [z(y)} = 
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Volume distance: 
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